A Numerical Renormalization Group for Continuum One-Dimensional Systems 



O 

o 

(N 
S3 

a 

(N 



Robert M. Konik and Yury Adamov 
CMPMS Department, Brookhaven National Laboratory, Upton, NY 11973 
(Dated: February 6, 2008) 

We present a renormalization group (RG) procedure which works naturally on a wide class of 
interacting one-dimension models based on perturbed (possibly strongly) continuum conformal and 
integrable models. This procedure integrates Kenneth Wilson's numerical renormalization group 
with Al. B. Zamolodchikov's truncated conformal spectrum approach. Key to the method is that 
such theories provide a set of completely understood eigenstates for which matrix elements can 
be exactly computed. In this procedure the RG flow of physical observables can be studied both 
numerically and analytically. To demonstrate the approach, we study the spectrum of a pair of 
coupled quantum Ising chains and correlation functions in a single quantum Ising chain in the 
presence of a magnetic field. 
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The numerical renormalization group (NRG), as de- 
veloped by Kenneth Wilson Q, is a tremendously suc- 
cessful technique for the study of generic quantum im- 
purity problems, systems where interactions are confined 
to a single point. But the NRG as such is not directly 
generalizable to systems where interactions are present 
in the bulk. The natural generalization of the NRG in 
real space treats boundary conditions between RG blocks 
inadequately, leading to qualitatively inaccurate results. 
To overcome this difficulty, Steven White developed the 
density matrix renormalization group (DMRG) [2] . This 
tool is now ubiquitous in the study of low dimensional 
strongly correlated lattice models and can access both 
static and dynamic quantities 

In this letter we offer a distinct generalization realiz- 
ing a renormalization group procedure for a wide range of 
strongly interacting continuum one-dimensional systems. 
It can treat any model which is representable as a confor- 
mal or integrable field theory (CFT/IFT), with Hamil- 
tonian, Hq , plus a relevant [J] perturbation (of arbitrary 
strength), H pclt . Beyond this there is no real constraint 
on Hq or H port . In particular, the full theory, Hq + H pclt , 
need not be integrable or conformal. Thus the technique 
can handle a standard array of models of perturbed Lut- 
tingcr liquids or Mott insulators. It also capable of treat- 
ing disordered systems, either by envisioning H pcrt as a 
random field or, equally well, considering a non-unitary 
supersymmetric CFT, Hq arising from disorder averag- 
ing a system with quenched disorder This technique 
also allows the study of coupled CFT/IFTs, allowing the 
study of systems between one and two dimensions. In all 
cases, the low energy spectrum and correlation functions 
of the model are computable. 

Our starting point is the truncated spectrum approach 
(TSA) pioneered by Al. B. Zamolodchikov. The TSA 
was developed to treat perturbations of simple confor- 
mal field theories. While straightforward in conception, 
it has an advantage over other numerical techniques in 
that it analytically embeds strongly correlated physics 
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FIG. 1: A schematic of the finite sized system, both in real 
space and in terms of energy levels, analyzed in the TSA 
procedure. 



at the start, dramatically lessening the computational 
burden. In one of the TSA's first applications, Al. B. 
Zamolodchikov studied a critical Ising chain in a mag- 
netic field [fl , a continuum version of the lattice model 
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where o~f are the standard Pauli matrices and i indexes 
the sites of the lattice. The continuum model, itself inte- 
grable, has a complicated spectrum of eight fundamental 
excitations. The TSA was able to produce the gaps of 
the first five excitations within 2% of the analytic, infi- 
nite volume values by diagonalizing a mere 39x39 matrix. 
To see how remarkable this is, consider that in a compu- 
tationally equivalent exact diagonalization of the lattice 
model, one would be limited to studying a five site chain. 

The TSA begins by taking the model to be studied, 
H = Hq + -ffpert, and placing it on a finite ring of cir- 
cumference, R. Doing so makes the spectrum discrete 
(see Fig. 1). We nonetheless expect to be able to ob- 
tain infinite volume results provided we work in a regime 
where RA 3> 1 with A a characteristic energy scale of 
the system. In the discrete system, the spectrum can 

then be ordered in energy, |2), Non-pcrturbative 

information is input in the next step of the TSA where 



the matrix elements, H p 



|-ffpert|i)j arc computed 
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exactly. It is important to stress this is always a practi- 
cal possibility. If Hq is a CFT, the attendant Virasoro 
algebra (or, equally good, some more involved algebraic 
structure such as a current or a W-algebra) permits the 
computation of such matrix elements straightforwardly. 
If Hq is instead an IFT, such matrix elements are avail- 
able through the form- factor bootstrap programme [HI ]. 
In an IFT, the matrix elements which we will want to 
focus on involve states, with few excitations and, as 
such, are readily computable. 

With the matrix elements in hand, one can then ex- 
press the full Hamiltonian as a matrix. The penultimate 
step in the procedure is to truncate the spectrum at some 
energy, £7trunc; making the matrix finite. This matrix 
is then numerically diagonalizcd from which the spec- 
trum and correlations functions can be extracted. When 
Hq represents a theory with a relatively simple set of 
cigenstates, this procedure, even with a crude truncation 
of states, works remarkable well in extracting the spec- 
trum. However when the starting point Hamiltonian, Hq, 
is more complicated (say a CFT based on a SU(2)k cur- 
rent algebra such as would be encountered in the study 
of spin chains), or one is interested in computing corre- 
lation functions in the full theory, Hq + H pclt , the simple 
truncation scheme ceases to produce accurate results at 
reasonable numerical cost. For example, errors in an ex- 
citation gap, A, introduced by the truncation of states 
behave as power laws, SA ~ -E'trunc ( a ~ !)■ Thus in- 
creasing E c does not dramatically reduce 6Q while at the 
same time greatly increasing the cost of the exact diag- 
onalization routine which should scale similarly to the 
partition of integers, i.e. as e l3 ^ Etr " nc /E tTnnc . It is the 
aim of this letter to outline an RG technique offering a 
dramatic improvement on this truncation scheme. 
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1) Take the first N+A N states of the theoryl 

2) Compute Hamiltonian and diagonalize 

3) Form a new basis of states using the first N eigenstates 
+ the next A N states in the original basis N NtVX 



4) Recompute Hamiltonian and numerically diagonalize 

5) Return to 3) 

FIG. 2: An outline of the NRG algorithm 

Our framework hews closely to the original Wilsonian 
conception of the NRG. In developing the NRG for the 
Kondo model, Wilson transformed the original Kondo 
Hamiltonian using the "Kondo basis" to a lattice model 
of an impurity situated at the end of a infinite half line 
with sites far from the end characterized by rapidly di- 
minishing matrix elements. We, in a sense, start in this 
position. The ordering in energy of the states provided 
by Hq is in direct analogy to the half-line on which the 
impurity lives in Wilson's Kondo work. The next step 
in Wilson's NRG is an iterative numerical procedure by 



which at each step a finite lattice is expanded by one 
site, the model diagonalized, and high energy eigenstates 
thrown away. It is this iterative procedure that wc mimic. 
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FIG. 3: Plots showing the behavior of A2 as a function of the 
truncation energy, Etrunc (f° r R = 5A _1 ). 

Let us denote the initial basis offered by Hq as 
We begin by keeping a certain number, say TV + ATV, 
of the lowest energy states, {|i)}^j[ AAr (in blue in Fig. 
2). We diagonalize the problem so extracting an ini- 
tial spectrum and set of eigenvectors, {\i)}f,Ji AN (in red 
in Fig 2). We then toss away a certain number, ATV, 
of the eigenvectors corresponding to the highest energy, 

1. e. {I^l^v+i ■ A new basis is then formed, consist- 
ing of the remaining eigenvectors together with the first 
ATV states of that we had previously ignored, 
he. {|i)}£Li U {IOlSv+AJV+1' an d tne procedure is re- 
peated. We present the technique schematically in Fig. 

2. Convergence of the procedure in the Kondo problem is 
promoted by the small matrix elements involving sites far 
from the impurity. While matrix elements in the proce- 
dure just described grow progressively smaller under the 
NRG (scaling as l/E tr unc), here numerical convergence is 
not necessarily the goal. Rather we aim to merely bring 
the quantity into a regime where its flow is governed by 
a simple flow equation. 

The algorithm just described implements a Wilsonian 
RG in reverse. It does so at all loop orders and so the 
RG flow it describes is exact. As the flow proceeds it, 
however, evolves closer and closer to a flow described by 
a one-loop equation. Analytically, the equation is nearly 
trivial as it is given solely in terms of the anomalous di- 
mension, Q!q, of the flowing quantity, Q (whether it be 
an energy eigenvalue or a matrix element). More specif- 
ically, 



dAQ 



d\nE t 



-g{a Q )AQ, 



(2) 



where AQ = Q(£Wnc) - Q(-Btrunc = 00) describes the 
deviation of the quantity as a function of the trunca- 
tion energy from its 'true' value (i.e. the value where 
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the cutoff in energy is taken to infinity). The function 
g{oLQ) can be determined exactiy using high energy per- 
turbation theory, well-controlled provided H pcrt is rele- 
vant. For example, if Q is some energy eigenvalue, then 
g(a Q ) = o>q = 1. 

The virtue of Eqn. @ is that it allows us to run the RG 
in two steps. We first implement the NRG as described 
above until we reach a truncation energy placing us safely 
in the one-loop regime. We then continue the RG by 
merely integrating the above equation allowing us to fully 
eliminate the effects of truncation. 

We now consider two examples using this RG proce- 
dure, one where we compute the spectrum of a model 
and one where we analyze correlation functions. Both 
examples are chosen so that a straight application of the 
TSA leads to poor results. 

Spectrum: In the first example, we consider a pair of 
quantum critical Ising chains coupled together: 

H = II,';-* 1 + H 1 ^ 2 ; H pelt = -\J dxatal (3) 

This model is known to be integrable and to have a spec- 
trum equivalent to the sine-Gordon model at /3 2 = tt 
[Io| . that is, a spectrum with a pair of solitons with gap, 
A s i = A S 2, together with a set of six bound states with 
gaps, Afe = 2A S sin(7rfc/14), k = 1, •••,6. By compar- 
ing conformal perturbation theory with a thermodynamic 
Bethe ansatz analysis, Agi^ can be expressed in terms 
of the coupling constant, A: A s i 2 = 11.2205920A with 
A = A 4 / 7 /(2^) 3 / 7 0. 

The underlying finite volume Hilbert space of Hq is 
considerably more complicated than that of single Ising 
chain. In a single chain there are four potential sectors of 
the Hilbert space a sector, /, composed of even num- 
bers of half-integer fcrmionic modes acting over a unique 
vacuum, \I); a sector, F, composed of odd numbers of 
half-integer modes over |/); and finally two sectors, cr/fi, 
composed of even numbers of both right and left moving 
integer fermionic modes over degenerate vacuua, |er)/|^}. 
The sectors a and [i are connected by applying a product 
of an odd number of even mode fermionic operators. Un- 
der periodic boundary conditions, however, the Hilbert 
space of a single chain is reduced to two sectors, / and a. 
In the two chain case, this is no longer true. Not only do 
we have sectors of the form I® I, I ®a, a®I : and a <E> 
(such tensor products arising naturally from considering 
two chains), but F®F, F®fi, n®F, and fi®(i. Unlike a 
single chain, all possible sectors are consistent with peri- 
odic boundary conditions. The Hilbert space that results 
is thus much larger and applying the TSA with a simple 
truncation scheme leads to poor results for the spectrum. 
Computing the spectrum of this model is thus an ideal 
testing ground for our proposed RG procedure. 

In Figure 3 we outline the procedure by which we 
extract the values of the spectrum of the two coupled 
Ising chains focusing for specificity on the second bound 
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A s i 


11.2206 


11.92/12.67 


11.32/11.54 


11.17(2)/11.15(5) 


A s2 


11.2206 


11.92/12.66 


11.32/11.54 


11.17(2)/11.16(6) 


Aj 


4.9936 


5.29/5.61 


5.03/5.12 


4.97(l)/4.97(2) 


A 2 


9.7369 


10.69/11.55 


9.89/10.24 


9.70(3)/9.7(l) 


A 3 


13.9918 


15.58/16.65 


14.33/14.84 


14.02(5)/14.20(5) 


A 4 


17.5452 


19.672/20.923 


18.69/18.03 


17.6(1)/17.7(1) 


A 5 


20.2188 


23.64/24.64 


20.80/21.62 


20.2(2) /20.5(2) 


A 6 


21.8785 


23.65/25.28 


22.39/23.08 


21.8(1)/21.8(2) 



TABLE I: The excitation energies for two coupled Ising chains 
at values of R = 4/5A -1 . (in units of A = A 4/7 /(2tt) 3/7 ). 



state, A2, in the spectrum. We first show the results 
of a straight TSA analysis (red squares) as a function 
of increasing truncation energies (given in units of 1/-R)- 
While the gap, A2, is converging towards its infinite vol- 
ume value, 9.7368648 ... A, it is doing so only slowly and 
at exponentially increasing numerical cost. We have per- 
formed the straight TSA analysis up to level 15 (i.e. keep- 
ing states with energies less than 15/ i?) where the Hilbert 
space contains r* 3500 states. At this point, the TSA pro- 
duces a result deviating by « 20% from the exact value. 
We also plot the value of A2 as given by the NRG algo- 
rithm as it iterates through states of ever higher energy. 
Here we have run the algorithm so as to take into ac- 
count states up to level 25 (in total « 150000/sector). 
We see, reassuringly, that where the TSA results exists, 
the NRG algorithm produces matching results (at a frac- 
tion of the numerical cost). The NRG algorithm ends up 
producing a value of A2 with a 5% error. Finally we plot 
the result of fitting Eqn. ([2]) to the NRG results between 
level 20 and 25 (where we believe a one-loop RG equa- 
tion describes the NRG flow). Extrapolating the fit to 
Etraac = 00 gives A2 = 9.70A, a value deviating from the 
exact result by 0.5%. 

In Table 1 we present the results of our RG analysis 
on the complete spectrum of the two coupled chains. In 
the first column we provide the exact value of spectrum 
as determined by integrability and the TBA analysis of 
Rcf. [7j. In the second column we give the values of 
the spectrum at two different system sizes (R = 4/5A -1 ) 
computed using a straight TSA analysis truncating at 
level 10 (i.e. keeping w 600 states in each of the relevant 
sectors). We see that the disagreement with the exact 
result ranges up to 20%. In the third column we give the 
results coming from applying the NRG algorithm (again 
iterating until we have reached level 25). We see a marked 
improvement over the TSA analysis, but still we obtain 
results with errors ranging up to 5%. In the final column, 
we give the results for the spectrum arrived at by fitting 
the one- loop RG equation (Eqii.®) to the NRG data. 
We see that our errors are now less than 1% . 
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FIG. 4: A plot of the matrix element, (0\a(0)\A 1 (p)A 1 (-p)) 
as a function of energy, uj = 2(p 2 + mf) 1 ^ 2 . 

Correlation Functions: We now turn to the compu- 
tation of correlation functions using the above described 
RG methodology. For simplicity we consider only T = 
response functions although a generalization to finite 
temperature multi-point functions is readily realizable. 
At zero temperature, the imaginary piece of a retarded 
correlation function, G rc t(x,t) = (0(x, £)0(O)) rc t, has a 
spectral decomposition, Soo(x,uj) = — ilmG rc t(x, uj > 
0) equal to [H 



Sooix, u) 



*K0|O(0,0)|n; S „)| 2 <5(cc,-£; s J,(4) 



where \n; s n ) is an eigenstate of the system with en- 
ergy/momentum E Sn /P Sn built out of n fundamental sin- 
gle particle excitations carrying internal quantum num- 
bers s n . Thus the computation of any response func- 
tion is equivalent to the computation of a number of ma- 
trix elements of the form (O|0(x, 0)|n; s n ). Ostensibly to 
compute the response function fully one would need to 
compute an infinite number of such matrix elements. In 
practice if one is interested in the response function at 
low energies only a small finite number of such matrix 
elements need be computed [l2j. We will illustrate the 
computation under the RG of one such non-trivial matrix 
element for a critical Ising chain in a magnetic field (i.e. 
Eqn. (flj ) . While the spectrum of this theory rapidly con- 
verges upon the increase of -Etruncj the matrix elements 
are less well behaved. And unlike the two-chain case, 
analytical results are available for comparison Q. Thus 
this computation is a good test of our RG methodology. 

We specifically compute the two excitation contribu- 
tion f 2 (oj) = (0|cr(0)|Ai(p)Ai(-p)|0) with p = (uj 2 - 
A 2 j), to the spin-spin correlator, S a a(p = 0,w > 0). 
(Here Au is the gap of the lowest lying single particle 
excitation, A\ , in an Ising chain in a magnetic field.) This 
contribution takes the form 



We are able to compute the necessary infinite volume 
matrix element over a continuous range of energies by 
studying a single matrix element in finite volume where 
the spectrum is discrete. We do so by continuously vary- 
ing the system size, R. Under such variations, the energy 
uj = 2(p 2 + nil) 1 / 2 of the state, \Ai{p)Ai(— p)), changes 
continuously due to the quantization condition of the mo- 
mentum, p, (i.e. p = 2im/R + S(p, —p) where 5(p, —p) is 
a two-body scattering phase). 

In Figure 4 by varying R we parametrically plot the 
results of our computations of f^ (p) vs its exact value [§] . 
A straight application of the TSA (with a level 10 trunca- 
tion) produces acceptable results at higher energies but 
does poorly at energies around threshold, 2Ai. At larger 
values of R (and so smaller energies), the TSA breaks 
down. The TSA curve in this region is then double val- 
ued. Computing the same matrix element with the NRG 
algorithm leads to a considerable improvement but at the 
lowest energies a deviation from the exact result remains 
(see inset to Fig. 4). RG improving the computation of 
/2(f) largely removes this discrepancy even at energies 
next to threshold. (In applying Eqn. (2) to Q = f2(p), 
perturbation theory yields, g(aj 2 ) = 2(1 — 1/8) where 
1/8 is the anomalous dimension of the spin operator, a). 

In conclusion we have presented an RG scheme by 
which a large number of one-dimensional continuum 
models can be studied with quantitative accuracy. With 
this methodology, both the spectrum and spectral func- 
tions of a model can be determined. 
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with A. Tsvelik. 
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